Lowest Order Constrained Variational calculation for Polarized 
Liquid ^He at Finite Temperature 

G.H. Bordbar* \ M.J. Karimi and J. Vahedi 

Department of Physics, Shiraz University, Shiraz 71454, Iran 

Abstract 

We have investigated some of the thermodynamic properties of spin polarized liquid ^He at 
finite temperature using the lowest order constrained variational method. For this system, the 
free energy, entropy and pressure are calculated for different values of the density, temperature 
and polarization. We have also presented the dependence of specific heat, saturation density and 
incompressibility on the temperature and polarization. 
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I. INTRODUCTION 



The spin polarized liquid ^He at finite temperature is an interesting many-body system 
which obeys the Fermi-Dirac statistics. Several theoretical techniques and experimental 
method have been used for investigating the properties of unpolarized liquid ^He 
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mm. 



Recently, we have studied the unpolarized hquid ^He at finite temperature and 
the spin polarized liquid '^He at zero temperature j^, Q, 11|. In this calculations, we have 
applied the lowest order constrained variational method based on the cluster expansion of the 



energy functional. This method is a power 



calculations of dense matter 
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ul microscopic technique used in the many-body 
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In the present work, we use the lowest order constrained variational method to investigate 
some of the thermodynamic properties of spin polarized liquid "^e at finite temperature. In 
our calculations, we employ the Lennard- Jones 19| and Aziz [20] potentials. 



II. METHOD 



We consider a system consisting of interacting ^He atoms with A^"*" spin up and spin 
down atoms. These atoms have a different chemical potential (/u) and different Fermi-Dirac 



distribution function 



23| as follows, 



n^ik) = , (1) 

exp[/3(^-^±)] + l' 

where P = and T is the temperature. We define the total number density by p and the 
spin asymmetry parameter hj ^ = is unpolarized and ^ = 1 is fully polarized cases), 

P = P^ + P', 

« = (2) 

Thus, to determine the chemical potential for every p, T and ^, we solve the following 
equation by numerical method, 

1 POD 

P^ = 7^ n^{k)kHk. (3) 

We calculate the energy using the lowest order constrained variational method based on the 
cluster expansion of the energy. We consider up to the two-body energy term in cluster 



expansion, 

E = El + E2. 

For this system, the one-body energy per particle, Ei, is given by 



(4) 



El = Et + E{ 



Armr'^p 

The two-body energy per particle, E2, is defined 

1 



00 ^00 
Jo 



(5) 



where 



E2 = ^j:i'j\Wil2)\zj-ji), 



Wiru) = -{yf{ri2)f + f\ri2)V{Ti2). 
m 



(6) 



(7) 



In the above equation, V{ri2) is the two body potential between helium atoms and /(ri2) is 
the two-body correlation function. By considering the \i) as a plane wave, we have derived 
the following relation for the two-body energy per particle. 



E2 = -p I df 



1 - 



Wir), (8) 



where 



sin(fcr) ±.,.,2j; 
-n^{k)k^dk. 



/o kr 

Now, we minimize the two-body energy Eq. ([H]) with respect to t 
body correlation function subject to the normalization constraint 



(9) 



le variations in the two- 
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Y.m\n2) - f\n2)\^J - ji) = 1. 



(10) 



The normalization constraint is conveniently re-written in the integral form as 



p / dr[h^{r) - P{r)]T{r) 



where 



Tir) 



and h{r) is the Pauli function, 



A-n^p^ 
h{r) = [r(r)] 



(11) 



(12) 



(13) 
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The minimization of the two-body Energy E2 gives the following Euler-Lagrange differential 
equation for the two-body correlation function 



/"(r) + /(r) 



r(r) 



m 



(A + y(r))/(r) = 0. 



(14) 



r(r) 

The lagrange multiplier A imposes by normalization constraint. The two-body correlation 
function calculated by numerically integrating Eq. ( fT4|) and then the energy per particle of 
system can be obtained. Finally the Helmholtz free energy per particle can be determined 
using the following equation, 

F = E-TS. (15) 



S is the entropy per particle 23|] which is given by 



where 



1 



s = s+ + s- 



n^{k) ln(n^(A;)) + (1 - n^{k) ln(l - n^{k)) 



k^dk. 



(16) 



(17) 



27r2p Jo 

By calculating E, S and F, we can obtain the thermodynamic properties of the spin polarized 
liquid ^He at finite temperature. 



III. RESULTS 



The two-body correlation function at T = l.OK and T = A.OK for different values of 
the spin asymmetry parameter {^) are shown in Fig. [H This figure shows that at low 
temperatures, the correlation function increases more rapidly and reaches the limiting value 
(/(r) = 1) at larger interatomic distances (r). It is shown that by increasing ^, the correlation 
function reaches the limiting value at lower r. It is seen that our results for the correlation 
function with the Lennard- Jones and Aziz potentials are identical. 

Our results for the free energy of liquid ^He with the Lennard- Jones and Aziz potentials 
at different spin asymmetry parameters {^) are presented in Fig. [2l We can see that at 
low densities (p < 0.011A~^), the free energies with the Lennard- Jones and Aziz potentials 
are similar. However, at higher densities, the difference between these free energies are 
noticeable. Fig. [2] shows that the free energy increases by increasing spin asymmetry 
parameter and this variation becomes higher at high temperatures. We see that the free 
energy doesn't show a bound state above a certain value of temperature (called flashing 
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temperature, Tf ). We have found Tf for different values of ^ which given in Table HI From 
Fig. [21 it is also seen that for each value of spin asymmetry parameter, the density of 
saturation point (minimum point of the free energy) decreases by increasing temperature. 
We see that the saturation density (po) increases by increasing spin asymmetry parameter 
for each value of temperature. Our results for the saturation density are given in Table ITTl 
Table [T] (Table [Tll) shows that for all values of ^ the flashing temperature (the saturation 
density) with the Aziz potential is lower than that of the Lennard- Jones potential. 

The entropy of liquid '^He is shown as a function of density for two different tempera- 
tures and various asymmetry parameters in Fig. [3l It is seen that the entropy decreases 
by increasing both density and spin asymmetry parameter and increases by increasing tem- 
perature. We have found that for each temperature, the difference between entropy of fully 
polarized and unpolarized liquid ^He decreases for higher values of density, especially at low 
temperatures. 

In Figs, mini we have shown the free energy, entropy and the specific heat of liquid ^He 
as a function of temperature for different values of at p = 0.0166^4^^. In these figures, 
the experimental results l|] for the unpolarized liquid ^He are also shown for comparison. 
Fig. m indicates a good agreement between our results for the free energy with the Lennard- 
Jones potential and the experimental results. From Figs. HHHl it is seen that the free energy 
(entropy and specific heat) decreases (increase) by increasing temperature. In addition, 
we see that for the free energy, entropy and the specific heat of liquid ^He, the differences 
between fully polarized and unpolarized cases increase by increasing temperature. 

The isothermal pressure is obtained from the free energy using the following equation. 



P(p,T)=p2— (18) 

In Fig. U\ we have plotted the isothermal pressure as a function of density for the various 
temperatures and ^. It is seen that the equation of state with the Aziz potential is stiffer than 
that of the Lennard- Jones potential. This figure shows that at low temperatures we have a 
liquid-gas phase equilibrium for the liquid ^He. We know that at the critical temperature (T^) 
there is no liquid-gas phase equilibrium. Our results of the critical temperature for different 
values of ^ are presented in Table IIIII We see that the critical temperature increases by 
increasing ^. For all ^ it is seen that the (Tc) with the Aziz potential is lower than that of 
the Leonard- Jones potential. 
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For liquid ^He, the saturation incompressibility at saturation density is given by, 



i^o(T) =9[p- 




) 



(19) 



In Table [TTl our calculated values of the saturation incompressibility are presented for differ- 
ent values of ^ and T. It is seen that for each value of ^, the incompressibility decreases by 
increasing temperature and at each temperature, it increases by increasing ^. From Table 
im we see that for all values of spin asymmetry parameter our results for the Kq with the 
Aziz potential are lower than those of the Leonard- Jones potential. 

IV. SUMMARY AND CONCLUSION 

We have considered a system consisting of the Helium atoms ('^He) with an asymmetrical 
spin configuration. For this system, we have calculated some of the thermodynamic proper- 
ties using the lowest order constrained variational method with the Lennard- Jones and Aziz 
potentials. Our calculations lead to the following conclusions for the liquid ^He: 

• At high temperatures, the correlation function reaches the limiting value for the smaller 
values of the interatomic distances. 

• The results of the correlation function with the Lennard- Jones and Aziz potentials are 
very similar. 

• The free energy with the the Lennard- Jones and Aziz potentials are identical for small 
values of the density. 

• The free energy, saturation density and incompressibility decrease by increasing tem- 
perature. 

• The free energy, saturation density, incompressibility and critical temperature increase 
by increasing spin asymmetry parameter. 

• The entropy, specific heat and flashing temperature decrease by increasing spin asym- 
metry parameter. 

• The entropy decreases by increasing density and increases by increasing temperature. 
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• The equation of state with the Aziz potential is stiffer than that of the Lennard- Jones 
potential. 

• For all values of spin asymmetry parameter our results for the flashing temperature, 
saturation density, critical temperature and saturation incompressibility with the Aziz 
potential are lower than those of the Leonard- Jones potential. 
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TABLE I: The flashing temperature (Tf) for different values of ^. 





Lennard- Jones 


Aziz 




TfiK) 


Tf{K) 


0.0 


2.8 


2.6 


0.4 


2.7 


2.5 


0.8 


2.6 


2.3 


1.0 


2.2 


1.9 



TABLE II: The saturation density (po) and the saturation incompressibihty (Kq) of hquid ^He for 
different values of ^ and temperature. 



Lenard- Jones Aziz 





T{K) 




Ko{K) 




Ko{K) 


0.0 


0.5 


0.0117 


185 


0.011 


174 


0.0 


1.0 


0.0116 


163 


0.0109 


157 


0.0 


1.5 


0.0113 


135 


0.0106 


121 


0.4 


0.5 


0.0117 


187 


0.011 


183 


0.4 


1.0 


0.0116 


172 


0.0109 


162 


0.4 


1.5 


0.0114 


145 


0.0107 


132 


0.8 


0.5 


0.0120 


237 


0.0113 


212 


0.8 


1.0 


0.0120 


232 


0.0113 


206 


0.8 


1.5 


0.0119 


203 


0.0112 


179 


1.0 


0.5 


0.0128 


275 


0.0119 


260 


1.0 


1.0 


0.0128 


263 


0.0119 


246 


1.0 


1.5 


0.0127 


212 


0.0118 


201 
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TABLE III: The critical temperature (Tg) for different values of ^. 





Lennard- Jones 


Aziz 




Tc{K) 




0.0 


5.3 


5.1 


0.4 


5.5 


5.3 


0.8 


6.2 


5.8 


1.0 


6.5 


6.1 
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1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



r (A) r (A) 



FIG. 1: The correlation function versus the interatomic distance (r) with the Leonard- Jones (full 
curves) and Aziz (dotted curves) potentials at p = 0.0166^-3 for ^ = 0.0 (a), ^ = 0.4 (b), ^ = 0.8 
(c) and i = 1.0 (d). 
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^ 1 1 1 1 1 1 1 1 1 1 1 1 

0.005 0.007 0.009 0,011 0.013 0.015 0.017 



Density (A ) 



-14 I • 1 • 1 • 1 • 1 • 1 • 1 

0.005 0.007 0.009 0.011 0.013 0.015 0.017 



Density (A'') 




0.005 0.007 0.009 0.011 0.013 0.015 0.017 0.005 0.007 0.009 0.011 0.013 0.015 0.017 

Density (A'' ) Density (A'' ) 



FIG. 2: The free energy of liquid ^He versus density at different values of temperature (T = 
1.0, 2.0, 3.0,4.0, 5. OiC)witli the Leonard-Jones (full curves) and Aziz (dotted curves) potentials for 
i = 0.0 (a), ^ = 0.4 (b), C = 0.8 (c) and ^ = 1.0 (d). 
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FIG. 4: The free energy of liquid ^He versus temperature for p = O.OIGGA"^ at different values 
of ^ with the Leonard-Jones (full curves) and Aziz (dotted curves) potentials. The experimental 
results (dashed curve) for unpolarized liquid '^He [l| are also shown for comparison. 
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p=0.0166A"' 


^=0.0 




yf ^=0-4 


- 


^=0.8 






[S^^^, 1,1, 


1,1, 



0.0 0.5 1.0 1.5 2.0 

Temperature (K) 



2.5 



FIG. 6: As Fig. HJ but for the specific heat. 
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FIG. 7: The pressure versus density at different values of temperature (T = 1.0, 3.0, 5.0, 6.0K) 
with the Leonard- Jones (full curves) and Aziz (dotted curves) potentials for ^ = 0.0 (a), ^ = 0.4 
(b), C = 0.8 (c) and ^ = 1.0 (d). 
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